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We analyze by exact Renormalization Group (RG) methods the infrared properties 
of an effective model of graphene, in which two-dimensional massless Dirac fermions 
propagating with a velocity smaller than the speed of light interact with a three- 
dimensional quantum electromagnetic field. The fermionic correlation functions are 
written as series in the running coupling constants, with finite coefficients that admit 
explicit bounds at all orders. The implementation of Ward Identities in the RG 
scheme implies that the effective charges tend to a line of fixed points. At small 
momenta, the quasi-particle weight tends to zero and the effective Fermi velocity 
tends to a finite value. These limits are approached with a power law behavior 
characterized by non-universal critical exponents. 

1. INTRODUCTION AND MAIN RESULT 

The charge carriers in graphene at half filling are effectively described by massless Dirac 
fermions constrained to move on a two-dimensional (2D) manifold embedded in three- 
dimensional (3D) space [9], with a Fermi velocity v that is approximately 300 times smaller 
than the speed of light. As a consequence, already without taking into account the interac- 
tions, the system displays highly unusual features as compared to standard 2D electron gases, 
such as an anomalous integer quantum Hall effect and the insensitivity to disordered- induced 
localization; most of these effects have already been experimentally observed [27, 28]. The 
study of many-body interactions among the charge carriers in graphene is of course very im- 
portant, particularly in view of recent experiments that suggest their relevant role in several 
physical properties of graphene [8, 22, 24, 34]. 

The effect of a weak short range interaction in graphene is quite well understood: it turns 
out that the behavior of the ground state is qualitatively similar to the free one, except 
that the Fermi velocity and the wave function renormalization are renormalized by a finite 
amount. This was expected on the basis of a power counting analysis [19-21]; recently, it 
has been rigorously proven in [16], where the convergence of the perturbative series was 
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established, by using the methods of constructive Quantum Field Theory (QFT) and by 
taking into full account the lattice effects (i.e., by considering the Hubbard model on the 
honeycomb lattice). 

The situation in the presence of long range interactions is much more subtle and still not 
completely understood. Their effect in graphcnc is usually studied in terms of a model of 
Dirac fermions interacting via a static Coulomb potential; retardation effects of the electro- 
magnetic (e.m.) field are neglected because the free Fermi velocity v is 300 times smaller 
than the speed of light c. At weak coupling, a logarithmic divergence of the effective Fermi 
velocity v{\s.) at the Fermi points and a finite quasi-particles weight have been predicted, 
on the basis of one-loop [18] and two-loops [26] computations. An unbounded growth of 
the effective Fermi velocity was also confirmed by an analysis based on a large- A'^ expansion 
[23, 33], which predicted a power law divergence of i'(k) at p^. It is not clear how to rec- 
oncile the logarithmic divergence expected from two-loops perturbative computations with 
the power law behavior found by large- A" expansions; moreover, the description in terms of 
Dirac fermions introduces spurious ultraviolet divergences that can produce ambiguities in 
the physical predictions [21, 26]. These difficulties may be related to a basic inadequacy of 
the effective model of Dirac fermions with static Coulomb interactions: the fact that the 
effective Fermi velocity diverges at the Fermi points (as predicted by all the analyses of 
the model) signals that its physical validity breaks down at the infrared scale where v{k.) 
becomes comparable with the speed of light; at lower scales, retardation effects must be 
taken into account, as first proposed in [17], where a model of massless Dirac fermions prop- 
agating with speed f ^ c and interacting with an e.m. field was considered. In [17] it was 
found that, at small momenta, the wave function renormalization diverges as a power law; 
this implies that the ground state correlations have an anomalous decay at large distances. 
Moreover, it was found that the interacting Fermi velocity increases up to the speed of light, 
again with an anomalous power law behavior. Despite its interest, the model proposed in 
[17] has not been considered further. The results in [17] were found on the basis of one-loop 
computations and in the presence of an ultraviolet dimensional regularization scheme. It is 
interesting to investigate whether the predictions of [17] remain valid even if higher orders 
corrections are taken into account and in the presence of different regularization schemes 
closer to the lattice cut-off that is truly present in actual graphene. 

The model we consider describes massless Dirac fermions in 2 + 1 dimensions propagating 
with velocity v < c, and interacting with a 3 + 1 dimensional photon field in the Feynman 
gauge. We will not be concerned with the instantaneous case (c — )■ oo); therefore, from 
now on, for notational simplicity, we shall fix units such that h = c = 1. The model is 
very similar to the one in [17], the main difference being the choice of the ultraviolet cut- 
off: rather than considering dimensional regularization, in order to mimic the presence of 
an underlying lattice, we explicitly introduce a (fixed) ultraviolet momentum cutoff both 
in the electronic and photonic propagators. The correlations can be computed in terms of 
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derivatives of the following euclidean functional integral: 

with, setting x = {xq,x) and x = {xi,X2) (repeated indexes are summed; greek and latin 
labels run respectively from to 2, 1 to 2), 

V{A,i;) := / [ej^,x^^,x - i^^^^,x^m,x] (1-2) 

J A 

B{J,(f)):= / rfx [j^,xJ"M,x + </'x'0x + ^xV'x] , 

where A is a three dimensional box of volume |A| = with periodic boundary conditions 
(playing the role of an infrared cutoff, to be eventually removed), the couplings e, z/^ are 
real and i^i = 1^2] the couplings z/^ arc counterterms to be fixed so that the photon mass is 
vanishing in the deep infrared. Moreover, ip^, ipy^ are 4-components Grassmann spinors, and 
the ii-th component j^,x of the current is defined as: 

Jo,x = ^V'xToV'x , J^^iv V'xTV'x , (1-3) 

where 7^ are euclidean gamma matrices, satisfying the anticommutation relations {7^, — 
—25/^^1,. The symbol P{dip) denotes a Grassmann integration with propagator 

9'^'\-) ■■= [ 7^e--- ^^"^" + ^"f \ o(k) . (1.4) 
^ ^ ' J (27r)3 + ^ ^ ' 

where (27r)^'^ J is a shorthand for |A|^^ X]k=27rn/L with n G I?, and Xo(k) = x(|k|) plays 
the role of a prefixed ultraviolet cutoff (here x(t) is a non increasing function from R"*" 
to [0, 1] such that x(t) = 1 if t < 1 and x(t) = if t > M > 1). Finally, A^^x are gaussian 
variables and P{dA) is a gaussian integration with propagator 

X<0)(^\ .- f ,»px Xo(p) _ /■ ^P^P3 ^ipx Xo(p) 



^ ^"^-y (27r)3^ 2|p| -y (27r)^^ p^+pi" ^'"'^ 

We perform an analysis based on the methods of constructive Renormalization Group (RG) 
for non relativistic fermions, introduced in [3, 12] (see [4, 25, 30, 31] for updated intro- 
ductions), which have already been proved effective in the study of several low-dimensional 
critical systems, such as one-dimensional (ID) interacting fermions [3, 5, 6], 2D critical Ising 
and vertex models [2, 15], the 2D Hubbard model on the square lattice at positive temper- 
atures [10], interacting fermions with asymmetric Fermi surfaces [11] and the 2D Hubbard 
model on the honeycomb lattice [16], just to mention a few. Compared to other RG ap- 
proaches, such as those in [29, 32], the advantage of the constructive methods we adopt is 
that they allow us to get a rigorous and complete treatment of the effects of the cut-offs and 
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a full control on the perturbative expansion via explicit bounds at all orders; quite remark- 
ably, in certain cases, such as the ones treated in [6, 10, 16], these methods even provide a 
way to prove the convergence of the resummed perturbation theory. 

By using these methods, we construct a renormalized expansion, allowing us to express 
the Schwinger functions, from which the physical observables can be computed, as scries in 
the effective couplings (the effective charges and the effective photon masses, also called in 
the following the running coupling constants) , with finite coefficients at all orders, admitting 
explicit A^! bounds (see Theorem 2.1 in Section 2 below). If the effective couplings remain 
small in the infrared, informations obtained from our expansion by lowest order truncations 
are reliable at weak coupling. The importance of having an expansion with finite coefficients 
should not be underestimated; the naive perturbative expansion in the fine structure con- 
stant is plagued by logarithmic infrared divergences and higher orders are more and more 
divergent. 

Of course the renormalized expansion is useful only as long as the running coupling 
constants are small. In fact, we do prove that they remain small for all infrared scales, by 
implementing Ward Identities (Wis) in the RG flow, using a technique developed in [6] for 
the rigorous analysis of Luttinger liquids in situations where bosonization cannot be applied 
(e.g., in the presence of an underlying lattice and/or of non-linear bands). The Wis that we 
use are based on an approximate local gauge invariance, the exact gauge symmetry being 
broken by the ultraviolet cut-off; its presence produces corrections to the "naive" (formal) 
Wis, which can be resummed and, again, explicitly bounded at all orders. The resulting 
modified Wis imply that the effective charges tend to a line of fixed points, exactly as in 
ID Luttinger liquids. We note that this is one of the very few examples in which Luttinger 
liquid behavior is found in dimensions higher than 1. 

Let us denote by (...)= lim|A|^oo (• • •)a expectation value with respect to the inter- 
action (1.2) in the infinite volume limit; our main result can be informally stated as follows 
(more rigorous statements will be found below). 

Main result. There exists a choice of such that, for k small, 



B(h.),u^,ri,fj,Veff are expressed by series in the effective couplings with finite coefficients 
that admit N\-bounds at all orders. Moreover: (i) the first non-trivial contribution to -B(k) 
is of second order in e; (ii) the first non-trivial contribution to z/^ is of second order in e 
and positive; (Hi) the first non-trivial contributions to ri,fi,Veff are, respectively: 





where: 



Z(k) ~ |k 







2e2 



'eff - 



1 - F{v) 




(1.8) 



V. 



5 



with 




^0 — arctan^o 



+ 



1 arctan^o" 



6:= 



(1.9) 



2v^ Co J 



V 



Note that the theory is not Lorentz invariant, because v ^ c; moreover, gauge symmetry 
is broken by the presence of the ultraviolet momentum cut-off. These two facts produce 
unusual features as compared to standard QFT models. In particular, the momentum cut- 
off produces correction terms in the Wis, which can be rigorously bounded at all orders; 
despite these corrections, one can still use the Wis to prove that the beta function for the 
effective charges is asymptotically vanishing at all orders, so that the model admits a line 
of (non-trivial) fixed points. 

The lack of gauge invariance due to the ultraviolet momentum cut-off makes it neces- 
sary (as in [7]) to introduce positive counterterms to keep the photon mass equal to zero. 
Similarly, it implies that the effective couplings with the temporal and spatial components 
of the gauge field are different and that the effective Fermi velocity Veff is not equal to 
the speed of light. However, it is possible to introduce in the bare interaction two different 
charges, Cq and ei, describing the couplings of the photon field with the temporal and spatial 
components of the current, which can be tuned so that the dressed charges are equal and 



A more realistic model for single layer graphene could be obtained by considering tight 
binding electrons hopping on the honeycomb lattice, whose lattice currents are coupled to 
a 3D photon field. A Renormalization Group analysis similar to the one in the present 
paper could be repeated for the lattice model, by extending the formalism in [16]. If the 
lattice model is chosen in such a way that lattice gauge invariance is preserved, we expect 
that its photon mass counterterms are exactly zero and that its effective Fermi velocity is 
equal to the speed of light. In any case (i.e., both in the presence or in the absence of 
lattice gauge invariance) , we expect the lattice model to have the same infrared asymptotic 
behavior of the continuum model considered here, provided that the bare parameters e^, 
of the continuum model are properly tuned. 

Finally, let us comment about the possibility of providing a full non perturbative con- 
struction of the ground state of the present model or, possibly, of a more realistic model of 
tight binding electrons hopping on the honeycomb lattice and interacting with e.m. forces. 
In this paper, we express the physical obervables in terms of series in the running coupling 
constants (with bounded coefficients at all orders) and we show that the running coupling 
constants remain close to their initial value, thanks to Ward Identities and cancellations 
in the beta function. Thus, the usual problem that so far prevented the non-perturbative 
construction of the ground state of systems of interacting fermions in d > 1 with convex 
symmetric Fermi surface (namely, the presence of a beta function driving the infrared flow 
for the effective couplings out of the weak coupling regime) is absent in the present case. 



Veff 



1. 
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Therefore, a full non-perturbative construction of the ground state of the present model 
appears to be feasible, by using determinant bounds for the fermionic sector and cluster 
expansion techniques for the bosonic sector. Of course, the construction is expected to be 
much more difficult than the one in [6] or [16], due to the simultaneous presence of bosons 
and fermions; if one succeeded in providing it, it would represent the first rigorous example 
of anomalous Luttinger-liquid behavior in more than one dimension. 

The paper is organized as follows: in Section 2 we describe how to evaluate the functional 
integrals defining the partition function and the correlations of our model in terms of an 
exact RG scheme (details are discussed in Appendix A and B); in Section 3 we describe 
the infrared flow of the effective couplings and prove the emergence of an effective Fermi 
velocity different from the speed of light (the explicit lowest order computations of the beta 
function are presented in Appendix C) ; in Section 4 we derive the Ward Identity allowing us 
to control the flow of the effective charges (details are discussed in Appendix D) and proving 
that the beta function for the charges is asymptotically vanishing; flnally, in Section 5 we 
draw the conclusions. 

2. RENORMALIZATION GROUP ANALYSIS 
A. The effective potential 

In this section we show how to evaluate the functional integral (1.1); the integration will 
be performed in an iterative way, starting from the momenta "close" to the ultraviolet cutoff 
moving towards smaller momentum scales. At the n-th step of the iteration the functional 
integral (1.1) is rewritten as an integral involving only the momenta smaller than a certain 
value, proportional to M~", where M > 1 is the same constant (to be chosen sufficiently 
close to 1) appearing in the deflnition of the cut-off function (see hues after (1.4)), and both 
the propagators and the interaction will be replaced by "effective" ones; they differ from 
their "bare" counterparts because the physical parameters appearing in their deflnitions (the 
Fermi velocity v, the charge e, and the "photon mass" i/^) are renormalized by the integration 
of the momenta on higher scales. In the following, it will be convenient to introduce the 
scale label h <0 as h :— —n. 

Setting X/i(k) := x(M~'*|k|), we start from the following identity: 



Xo(k) = Yl ' ^(k) X.(k) - X.-i(k) ; (2.1) 

h=—oo 

let i/j = E°=-oo V'^''^ and A = EL-oo^^''^ ^^^^^ {'ip^''^}h<o, {A'^^^}h<o are independent free 
fields with the same support of the functions fh introduced above. 

We evaluate the functional integral (1.1) by integrating the fields in an iterative way 
starting from A^°^; for simplicity, we start by treating the case J — (j) — 0. We 
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define V*^°^(A, ^) := V{A,tjj) and we want to inductively prove that after the integration of 
ip^^\A^'^\ -^C^+i), ^iC^+i) we can rewrite: 

gW(o,o) ^ g|A|E, J p(^^(</^))p(^^(<'^))eV"'UA(^''),yz^^(^'')) ^ (2.2) 
where P[dip^-^^) and P((i/1*^-^)) have propagators 

^ ' Z,(k) A;2 + i;,(k)2|^|2 ' 2|p| ' 

V^'^-* has the form 

2n ,-, m , n 



P2i 



n,m>0 p,fj, 
n+m>l 



m / m 2n \ 

•n^/^*-p^<U({ka'{Pi})^ EP^+E(-l)^k, , (2.4) 

i=l \j=l i=l / 

and Eh, ■^/i(k), ^^^(k) and the kernels Wmji^p^f^ will be defined recursively. 

In order to inductively prove (2.2), we split V^'') as CV^''^ + nV^^\ where 7^ = 1 - £ and 
jC, the localization operator, is a linear operator on functions of the form (2.4), defined by 
its action on the kernels Wm}n,p,n in the following way: 

>C<4(k) := <4(0) + kak<4(0) , (2.5) 
>C<4,(p,k) :=<4,(0,0), 

(P) := W^it (0) + pdMZiO) , riyij^pi, P2) := (0, 0) , 

and jCW^^ :— otherwise. As it will be clear from the dimensional analysis performed 
in Section 2 D below, these are the only terms that need renormalization; in particular, 
CWq'^i p{\s.) will contribute to the wave function renormalization and to the effective Fermi 

velocity, >C1^2 o^,/^(p) to the effective photon mass, and jCWi'^i p^^{p, k) to the effective charge. 
As a consequence of the symmetries of our model, see Appendix A, it turns out that 

<oUo) = o, i^S.(o,o) = o, <4(o) = o, 

W^ioUo) = -S,„,,M\,,, , apW^ifoUo) = (2.6) 



and, moreover, that 



V'k kak<? (O)V'k = -iz^,hk^tljy,l^.^l^i. (2.7) 
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with z^^h, ^n,h real, and zi^h — Z2,h, ^i,h — ^2,h- We can renormalize P{d'ip^-^^) by adding 
to the exponent of its gaussian weight the local part of the quadratic terms in the fermionic 
fields; we get that 

J P(rfV;(^^))P(dA(^^'^))e^''''(^'^^) = (2.8) 

where th takes into account the different normalization of the two functional integrals, V^'*^ 
is given by 

~ /■ f/k — 

V^''\A,^) = V^''\A,^) + J -^^iz,,nKi^^l,i^^ 

=: V('^H^,V')->Cv,V('^)(AV') , (2.9) 
and P{dip^-'^^) has propagator equal to 

^(</^)^k^ ^ Xfe(k) noko + ivh-i0<)k-^ 

with 

Zh-iO^) = ZhO^) + ZhZo,hXh(^) , (2.11) 
Zh-iik)vh-i{k) = Zh{k)vh{k) + ZhZi,hXhO^) ■ 

After this, defining := Z/j_i(0), we rescale the fermionic field so that 

v^^\A, ^ni^) = vW(A v^V') ; (2.12) 

therefore, setting 

Vh-\ := Vft-i(O) , eo,/i := ^^\h , ei,/jW,j_i = e2,ftW/,_i := -^^Xi^h , (2.13) 
we have that: 

£V('^)(A(^'^), ^ri^^^'^) = dx (Z,_,e„4|'^)4f ) - M'^.,,,4f )) , (2.14) 
where 

i^'^'W^'^ , il^'^^ ^^.-1 . (2.15) 

After this rescaling, we can rewrite (2.8) as 

J P(#(^'^))P(dA(^'^))e^''''(^'^'^) = el^l*" y" P(#(^'^-i))P(dA(^'^-^) • 

• J P(dV^W)P(d^W)e^"'n^'^''->+A(''),^(^(^''-)+^(''))) ^ (2.16) 
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where ijj'^-'^ ^\A^-^ have propagators given by (2.3) (with h replaced by /i — 1) and 
have propagators given by 

gW{k) ^ /fe(k)^7o/co + Wft_l(k)fc■7 ^(^^ ^ fh{p) 
Zh-l Zh-l fc2+^,.-i(k)2|^|2 ' 2|p| ' 

Mk) = xh{k) - xk-i{k) , f\{k) = /^A(k) . (2.17) 

At this point, we can integrate the scale h and, defining 

^v(^-^HA^z;rif)+WEh J p((^^W)p((^^W)gV('')(yi+A(''),^i;:i7(V.+^(''))) ^ ^2.18) 

our inductive assumption (2.2) is reproduced at the scale h — 1 with E^-i := Eh + th + Eh. 
Notice that (2.18) can be seen as a recursion relation for the efi^ective potential, since from 
(2.9), (2.12) it follows that 

V^^\A, ^/Zi^i^) = V^''\A, y/z^^/j) = V^'^XA, ^/Z~n^) - jC^V^'^^A, . (2.19) 

The integration in (2.18) is performed by expanding in scries the exponential in the r.h.s. 
(which involves interactions of any order in ip and A, as apparent from (2.4)), and integrating 
term by term with respect to the gaussian integration P{dip^^^)P{dA^'^^). This procedure 
gives rise to an expansion for the effective potentials V*-^-* (and to an analogous expansion 
for the correlations) in terms of the renormalized parameters {ef^_k,^fi,k, Zk-i,Vk-i}h<k<o, 
which can be conveniently represented as a sum over Feynman graphs according to rules 
that will be explained below. We will call {e^ ^, z/^ fc}/i<fc<o effective couplings or running 
coupling constants while {en^k}h<k<o are the effective charges 

Note that such renormalized expansion is significantly different from the power series 
expansion in the bare couplings e, i/^; while the latter is plagued by logarithmic divergences, 
the former is order by order finite. 

By comparing (1.1) and (1.2) with (2.2), (2.4) and (2.14), wc see that the integration of 
the fields living on momentum scales > M'^ produces an effective theory very similar to the 
original one, modulo the presence of a new propagator, involving a renormalized velocity 
Vh and a renormalized wave function Z^, and the presence of a modified interaction V*-^-*. 
The lack of Lorentz symmetry in our model (implied by the fact that v ^ 1) has two main 
effects: (1) the Fermi velocity has a non trivial flow; (2) the marginal terms in the effective 
potential are defined in terms of two charges, namely eo,/i and ei^h = G2,h, which are different, 
in general. 

B. Tree expansion 

The iterative integration procedure described above leads to a representation of the ef- 
fective potentials in terms of a sum over connected Feynman diagrams, as explained in the 
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following. The key formula, which we start from, is (2.18), which can be rewritten as 

= E l^^l {V^'\A^^'\ ^/Z^.^^^'hn) , (2.20) 

n>l 

with 8^ the truncated expectation on scale /i, defined as 

£l{X{A^'^\i^^^^);n) ^ log j P(ciV^'V(^^^'^)e'''^^'''''^*"^ (2.21) 

If X is graphically represented as a vertex with external lines A^^^ and the truncated 
expectation (2.21) can be represented as the sum over the Feynman diagrams obtained 
by contracting in all possible connected ways the lines exiting from n vertices of type X. 
Every contraction corresponds to a propagator on scale /i, as defined in (2.17). Since 
is related to V*^'^^ by a rescaling and a subtraction, see (2.9) and (2.12), Eq.(2.20) can 
be iterated until scale 0, and V^'^~^^ can be written as a sum over connected Feynman 
diagrams with lines on all possible scales between h and 0. The iteration of (2.20) induces a 
natural hierarchical organization of the scale labels of every Feynman diagram, which will be 
conveniently represented in terms of tree diagrams. In fact, let us rewrite V^'*^ in the r.h.s. 
of (2.20) as VW(A, y/Z^xil)) = ZV^^\A, ^/Z^ijj) +nV^''\A, VZ^i/j), where Z -.^jC-jC^, see 
(2.9). Let us graphically represent V^'^\ CV^'^^ and TZV^'^^ as in the first hne of Fig.l, and let 
us represent Eq.(2.20) as in the second line of Fig.l; in the second line, the node on scale h 
represents the action of S^. Iterating the graphical equation in Fig.l up to scale 0, we end 

V(h) = « , Zv*> = . , 31 V*' = ^ 



V^i) = . . + . ^ 





h-l h h-1 h h-1 h \ h-1 



FIG. 1: Graphical interpretation of Eq.(2.20). The graphical equations for CV^'^~^\ TZV^^~^^ are 
obtained from the equation in the second line by putting an C, TZ label, respectively, over the 
vertices on scale h. 

up with a representation of V^'^^ in terms of a sum over Gallavotti-Nicold trees r [3, 13, 14]: 

V^''\A^^''\y/%7p^^''^) = Y^ V^^\t) , (2.22) 



N>iT&rh 



h,N 



where Th,N is the set of rooted trees with root r on scale hj. — h and N endpoints, see 
Fig. 2. The tree value V^^\t) can be evaluated in terms of a sum over connected Feynman 
diagrams, defined by the following rules. 
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FIG. 2: The effective potential V^'^^^^ can be represented as a sum over Gallavotti - Nicold trees. 
Tfie black dots will be called vertices of the tree. All the vertices except the first {i.e. the one on 
scale h) have an TZ label attached, which means that they correspond to the action of Tl£^^ , while 
the first represents The endpoints correspond to the graph elements in Fig.3 associated to the 
two terms in (2.14). 




V) 



FIG. 3: The two possible graph elements associated to the endpoints of a tree, corresponding to 
the two terms in the r.h.s. of (2.14). 



With each endpoint w of r we associate a graph element of type e or v, corresponding to 
the two terms in the r.h.s. of (2.14), see Fig. 3. We introduce a field label f to distinguish 
the fields associated to the graph elements e and u (any field label can be either of type 
A or of type -0); the set of field labels associated with the endpoint v will be called 1^. 
Analogously, if v is not an endpoint, we call the set of field labels associated with the 
endpoints following the vertex v on t. 

We start by looking at the graph elements corresponding to endpoints of scale 1: we group 
them in clusters, each cluster being the set of endpoints attached to the same vertex v 
of scale 0, to be graphically represented by a box enclosing its elements. For any Gu of scale 
(associated to a vertex v of scale that is not an endpoint), we contract in pairs some 
of the fields in Uw(=GvIw, in such a way that after the contraction the elements of G^ are 
connected; each contraction produces a propagator g^^^ or w'^'^\ depending on whether the 
two fields are of type ip or of type A. We denote by the set of contracted fields inside the 
box G^ and by = 1^, the set of external fields of G^; if v is not the vertex immediatly 
following the root we attach a label TZ over the box G^, which means that the TZ operator. 
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defined after (2.4), acts on the value of the graph contained in Gy. Next, we group together 
the scale-0 clusters into scale-(-l) clusters, each scale-(-l) cluster being a set of scale-0 
clusters attached to the same vertex v of scale —1, to be graphically represented by a box 
enclosing its elements, see Fig.4. 




-2-1 1 



FIG. 4: A possible Feynman diagram contributing to and its cluster structure. 

Again, for each v of scale —1 that is not an endpoint, if we denote by Wi, . . . , t'.^,, the 
vertices immediately following v on r, we contract some of the fields of Uj^^^^P^. in pairs, in 
such a way that after the contraction the boxes associated to the scale-0 clusters contained 
in G^ are connected; each contraction produces a propagator g'^~^^ or w^~^\ We denote by 
X„ the set of fields in U|Z,iP„. contracted at this second step and by — U|Z,iP„. \ the 
set of fields external to G^; iiv is not the vertex immediatly following the root we attach a 
label TZ over the box G^. 

Now, we iterate the construction, producing a sequence of boxes into boxes, hierarchi- 
cally arranged with the same partial ordering as the tree r. Each box G^ is associated to 
many different Feynman (sub-) diagrams, constructed by contracting in pairs some of the 
lines external to G^., with Vj, i = 1, . . . , s^,, the vertices immediately following v on r; the 
contractions are made in such a way that the clusters G.^^ , ■ ■ ■ , G^^^ are connected through 
propagators of scale hy. We denote by and by P^ the set of fields of type A and ip, 
respectively, external to Gy. The set of connected Feynman diagrams compatible with this 
hierarchical cluster structure will be denoted by r(r). Given these definitions, we can write: 

v«M-E/n^n^Va,(c;). 

Val(^) = [ n ^g/'ip,] [ n \^^g!;(/)]^^(^o)Val(^) , (2.23) 

m(g)=(-irf n (t-;)^f^[(n«5'"')( n 

V not e.p. '^■"'^ 



& V e.p. 

v*>v, 
hy*=hv+l 



where: (— l)'^ is the sign of the permutation necessary to bring the contracted fermionic fields 
next to each other; in the product over / e P^, can be either ip or ip, depending on the 
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specific field label /; S{vo) — ^[^ Z^/eP-^ P/ ~ S/eP''' i~^y^'^^^fj ' where e{f) — ± depending 
on whether ip is equal to ^ or ^; the integral in the third line runs over the independent 
loop momenta; s„ is the number of vertices immediately following v on t; TZ — 1 — C is the 
operator defined in (2.5) and preceding lines); — ii v — Vq, and otherwise — 1; gf^^ is 
equal to g^'^^ or to ■u;*^'^) depending on the fermionic or bosonic nature of the line i, and i e v 

(k) 

means that £ is contained in the box Gy but not in any other smaller box; finally, K^J is 
the matrix associated to the endpoints v* on scale A; + 1 (given by ieo,fc7o if v* is of type (a) 
with label p = 0, by iej^kVk'jj if v* is of type e with label p = j E {1, 2}, or by — M^f^ ^ if 
V* is of type u. In (2.23) it is understood that the operators TZ act in the order induced by 
the tree ordering (i.e., starting from the endpoints and moving toward the root); moreover, 
the matrix structure of gf'^ is neglected, for simplicity of notations. 



C. An example of Feynman graph 

To be concrete, let us apply the rules described above in the evaluation of a simple 
Feynman graph Q arising in the tree expansion of V^'^"^-'. Let Q be the diagram in Fig. 5, 
associated to the tree r drawn in the left part of the figure; let us assume that the sets P„ 
of the external hues associated to the vertices of r are all assigned. Setting 




FIG. 5: A possible Feynman diagram contributing to V^'^ and its cluster structure. 



(2.24) 



we can write: 



Val(e?) = 



1 Zh Zh-l_2 -2 Tijh -J 



4!2! Zh-i Zn-2 



dp 



■n 



(27r)3 
(2^ 



«'^'Hp)l'7Miff^'Hk + p)- (2.25) 



where 7^ [F(k + p)] = F(k + p) - F(0) - (k + p) • VF(0) = lik^+ p^){K + Pu)d^d,F{\^*). 
Notice that the same Feynman graph appears in the evaluation of other trees, which are 
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topologically equivalent to the one represented in the left part of Fig. 5 and that can be 
obtained from it by: (i) relabeling the fields in P^^^, P^^, (ii) relabeling the endpoints of the 
tree, (iii) exchanging the relative positions of the topologically different subtrees with root 
Vq. If one sums over all these trees, the resulting value one obtains is the one in Eq.(2.25) 
times a combinatorial factor 2^ ■ 3 ■ 4 (2^ is the number of ways for choosing the fields in 
and in P^^^; 3 is the number of ways in which one can associate the label u to one of the 
endpoints of scale /i + 1; 4 is the number of distinct unlabelled trees that can be obtained 
by exchanging the positions of the subtrees with root Vq). 

D. Dimensional bounds 

We are now ready to derive a general bound for the Feynman graphs produced by the 
multiscale integration. Let Wm,n'^p,^ be the contribution from trees with N end-points to the 
kernel Wm]n,p,n in 2.4, that is 

oo * oo 

<U({k.}, {P,}) = E E E ^K^) = E {p.}) , (2.26) 

N=iTeTh,N eer{T) iv=i 

l<l=2n 

where the * on the sum indicates the constraints that: U/gpA{py} = U'JLi{pj}; 

The iV-th order contribution to the kernel of the effective potential admits the following 
bound. 

Theorem 2.1 (N\ bound) Letih — maxft<fe<o{|e/i,fc|, |^;u,fe|} be small enough. IfZ^/Z^^i < 
e^^h and < v^-i < 1, for all h < k < and a suitable constant C > 0, then 

I KS„I I < (const.)^ef (f ) ! mM3— (2.27) 

where \\W^:i%\ \ := sup{k.},{p.} \ 

The factor 3 — 2n — m in (2.27) is referred to as the scaling dimension of the kernel with 2n 
external fermionic fields and m external bosonic fields; according to the usual RG thermi- 
nology, kernels with positive, vanishing or negative scaling dimensions are called relevant, 
marginal or irrelevant opeiaiois, respectively. Notice that, if we tried to expand the effective 
potential in terms of the bare couplings e, f^, the A^-th order contributions in this "naive" 
perturbation series could not be bounded uniformly in the scale h asm (2.27), but rather by 
the r.h.s. of (2.27) times \h\^ , an estimate which blows up order by order as h ^ — oo. 
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Proof. Using the bounds 

\\g^''\k)\\ < const • M-'^ , J dk\\g^^\k)\\ < const • M^'' , 

\w^^\k) \ < const • M-'' , J dk\w^^'\k) \ < const • M"^^ , (2.28) 

and the assumptions on Vk-i and Z^/Z^^i into (2.23), we find that, if r e 7^^^ and Q e r(T), 



|Val(^;)| < (const.)^e^- (2.29) 

V not e.p. V not e.p. 

v>vq 

where: n° is the number of propagators i E v, i.e., of propagators i contained in the box 
but not in any smaller cluster; is the number of vertices immediately following v on 
r; is the number of end-points of type u immediately following v on r (i.e., contained 
in but not in any smaller cluster); v' is the vertex immediately preceding f on r and 
= 2 if \P^\ = \P^\ = 2, z., = 1 is \P^^\ = 2\P^\ = 2 and z., = otherwise. The last 
product in (2.29) is due to the action of 71 on the vertices v > Vq that are not end-points. 
In fact, the operator TZ, when acting on a kernel W^'^{\p, k) associated to a vertex v with 
\P^\ = 2\P^\ = 2, extracts from wi\''^ the rest of first order in its Taylor expansion around 
p = k = 0: if \wi^''\p,k)\ < C{v), then |w£''\p,k)| = ^Kp^p + kak)M^g''^(p*, k*)| < 
(const. )M~'*''+'*'''C(v), where M"'*" is a bound for the derivative with respect to momenta 
on scale hy and M'^^' is a bound for the external momenta p, k; i.e., TZ is dimensionally 
equivalent to M~^^'"~^^'\ The same is true if TZ acts on a kernel W'3,o''^(pi, P2). Similarly, if 
TZ acts on a terms with |P^| = 2, it extracts the rest of second order in the Taylor expansion 
around k = 0, and it is dimensionally equivalent to k^d"^ ~ j\^-2(/i«-/i„/) ^s a result, we get 
(2.29). 

Now, let (n^) be the number of vertices of type e (of type v) following v on r. If we 
plug in (2.29) the identities 

^ (/i„-/i)(s,-l)= {hy-K){nl + n^,-l) 

V not e.p. V not e.p. 

V not e.p. V not e.p. 

Y {hv-h)m:= J2 (hv-K^K (2.30) 

V not e.p. V not e.p. 

we get the bound 

|Val(^?)|<(const.)^£t^M'^(^-l^''ol) £l!^^(h„-V)(3-|P.|-^.) . (2.31) 

^"0- V not e.p. 

v>vo 
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In the latter equation, 3 — |P^| is the scaling dimension of the cluster Gy, and 3 — \Py\ — Zy 
is its renormalized scaling dimension. Notice that the renormalization operator TZ has been 
introduced precisely to guarantee that ?> — \Pv\ — z,^ < for all by construction. This fact 
allows us to sum over the scale labels h < < 1, and to conclude that the perturbative 
expansion is well defined at any order N of the renormalized expansion. More precisely, 
the fact that the renormalized scaling dimensions are all negative implies, via a standard 
argument (see, e.g., [3, 14]), the following bound, valid for a suitable constant C (see (2.27) 
for a definition of the norm || • ||): 



from which, after counting the number of Feynman graphs contributing to the sum in (2.32), 
(2.27) follows. ■ 

An immediate corollary of the proof leading to (2.27) is that contributions from trees 
r G Th,N with a vertex v on scale = k > h admit an improved bound with respect to 
(2.27), of the form < (const. )^ej^(iV/2)! M'*(3-I-P"oI)m^(''-^), for any < ^ < 1; the factor 
jyj-${h-k) ^g^^ thought of as a dimensional gain with respect to the "basic" dimensional 
bound in (2.27). This improved bound is usually referred to as the short memory property 
(i.e., long trees are exponentially suppressed); it is due to the fact that the renormalized 
scaling dimensions — 3 — \Py\ — z^ in (2.31) are all negative, and can be obtained by 
taking a fraction of the factors M^^'"~^-"''>^'" associated to the branches of the tree r on the 
path connecting the vertex on scale k to the one on scale h. 

Remcirk. All the analysis above is based on the fact that the scaling dimension 3 — |P^| in 
(2.31) is independent of the number of endpoints of the tree r; i.e., the model is renormaliz- 
able. A rather different situation is found in the case of instantaneous Coulomb interactions, 
in which case the bosonic propagator is given by (2|p|)~^ rather than by (2|p|)^^. In this 
case, choosing the bosonic single scale propagator as w^'^\p) = Xoip)fhip){'2\p\)^^, one finds 
that the last bound in (2.28) is replaced by J dp < (const.) (dimensionally, this 

bound has a factor missing). Repeating the steps leading to (2.31), one finds a general 
bound valid at all orders, in which the new scaling dimension is 3 — \Py\ + + n^; this 
(pessimistic) general bound assumes that at each scale the loop lines of the graph are all 
bosonic. Perhaps, this bound can be improved, by taking into account the explicit structure 
of the expansion; however, it shows that the renormalizability of the instantaneous case, if 
true, does not follow from purely dimensional considerations and its proof will require the 
implementation of suitable cancellations. 




(2.32) 




j^{hv-h^,){3-\Pv\-Zv) 
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Syml")ol 


Doscripticni 


r 


Gallavotti - Nicolo (GN) tree. 


r 


Root label of the tree. 


Vq 


First vertex of the tree, immediately following the root. 


K 


Scale label of the tree vertex v. 


7h,N 


Set of GN trees with root on scale hr = h and with N endpoints. 


Iv 


Set of field labels associated with the endpoint of the tree v. 




Cluster associated with the tree vertex v. 


T 

-^v 


Set of contracted fields inside the box corresponding to the cluster . 


Pv 


Set of external fields of G^. 


Vi 


i-th vertex immediately following v on the tree. 


Sy 


Number of vertices immediately following the vertex v on the tree. 


P* 


Set of fields of type = A^tp external to G^- 


r(r) 


Set of connected Feynman diagrams compatible with the hierarchical cluster 




structure of the tree r. 




Number of propagators contained in G^ but not in any smaller cluster. 


< 


Number of end-points of type v immediately following v on the tree. 


v' 


Vertex immediately preceding v on the tree. 


nt 


Number of vertices of type jj^ = e,u following v on the tree. 




Improvement on the scaling dimension due to the renormalization. 



TABLE I: List of the symbols introduced in Sections 2B, 2D. 

E. The Schwinger functions 

A similar analysis can be performed for the 2-point function, see Appendix B. It turns 
out that, similarly to what we found above for the effective potentials, the 2-point function 
can be written in terms of a renormalized perturbative expansion in the effective couplings 
{e//,fe, T^fj,,k\k<o ^-nd in the renormalization constants {Z^, Vk}k<Oi with coefficients represented 
as sums of Feynman graphs, uniformly bounded as |A| oo; in contrast, the graphs forming 
the naive expansion in e, z/^j are plagued by logarithmic infrared divergences. 

More exphcitly, if M^' < |k| < M^+^, we get (see Eqs.(B.10)~(B.13)): 




(2.33) 



where S(k) is given by a formal power series in {e^^^, i^^,jt}fe<o with coefficients depending 
on {^jfc, i'fc}fe<0) s-nd starting from second order; under the same hypothesis of Theorem 2.1, 
the A'"-th order contribution to -B(k) is bounded by (const. )^(£_oo)'^(-?V/2)! uniformly in k. 
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Eq.(2.33) is equivalent to Eq.(1.7) of the main result (see the end of Section 3C below for 
the explicit relation between Zh,Vh and Z(k),v(k)). 

To prove our main result we need to control the flow of the effective charges at all 
orders in perturbation theory, and to do this we shall use Ward Identities, see Section 4. 
These are nontrivial relations for the three point functions, which can be related to the 
renormalized charges in the following way. Consider a theory with a bosonic infrared cutoff 
M^* , that is assume that the bare bosonic propagator is given by (1.5) with Xo{p) replaced by 
X[/i*,o](p) '■= Xo(p) ~ Xo{M~^*p), which is vanishing for |p| < M^* and it is equal to Xo(p) 
for |p| > M^*'^^; denote by (. . .)^. the expectation value in the presence of the bosonic 
infrared cutoff. As shown in Appendix B, setting eo,h '■= eo,/i, ei^h = (^2,h '■= Vh-iCi^h, and 
taking |q| = M^*, |q + p| < M^*, |p| <C M^* (we will be interested in the limit p — >■ 0), the 
following result holds (see Eqs.(B.14)-(B.16)): 

Om,P; V'q+pV'q)^* = ^^^*-l^(V'q+pV'q+p)^* + B^,,h*{p, Q)) (V'qV'q)^* , (2-34) 

where B^^h* is given by a formal power series in {efj_^k,i^fj.,k}h*<k<o, starting from second 
order and with the A^-th order of the series admitting a bound proportional to {eh*)^ {N/2)\, 
uniformly in k. Eq.(2.34) is one of the two desired equations relating the 3-point function to 
the 2-point function and the effective charge e^^h* . A second independent equation expressing 
the 3-point function in terms of the 2-point function and of the bare charge e will be derived 
in Section 4, see (4.4), using the (approximate) gauge invariance of the theory. Combining 
the two equations we will be able to relate e^,/i* to the bare charge e, for all h* < 0, and 
this will allows us to control the flow of the effective couphngs on all infrared scales. This 
procedure will be described in detail in the next two sections. 

3. THE FLOW OF THE EFFECTIVE COUPLINGS 

A. The Beta function 

A crucial point for the consistency of our approach is that the running coupling constants 
Gfj.,h,^fi,h are small for all h < 0, that the ratios Zh/Zh_i are close to 1, and the effective 
Fermi velocity Vh does not approach zero. Even if we do not prove the convergence of 
the series but only A^! bounds, we expect that our series gives meaningful information 
only as long as the running coupling constants satisfy these conditions. In this section we 
describe how to control their flow. We shall proceed by induction: we will flrst assume 
that £ = maxfe<o{|e^^fc|} is small, that Z^jZ^-x < e^^^ and < ^^/i < 1 for all /i < 
and a suitable constant C > 0, and we will show that, by properly choosing the values of 
the counterterms in (1.2), the constants Vf^^h remain small: max^<o{|t'^,/i|} < (const.) 
Next, once that the flow of Vfj^^h is controlled, we will study the flow of Z^ and Vh under the 
assumption that the constants e^^/j remain bounded and small for all /i < 0; we will show 
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that, asymptotically as /i — >■ — oo, Zy^ ~ M~^^, with rj — O(e^) a positive exponent, while Vh 
grows, approaching a limiting value fg// close to the speed of light. Finally, we shall start 
to discuss the remarkable cancellations following from a Ward Identity that guarantee that 
the constants e^^h remain bounded and small for all h < 0; the full proof of this fact will be 
postponed to Section 4 and Appendix D. 

The renormalized parameters obey to recursive equations induced by the previous con- 
struction; i.e., (2.6), (2.7), (2.11), (2.13) imply the flow equations: 

%i = l + Zo,/.:=l + ^L vn_^ = ^{vn + z,,h)■.^VH + Pl (3.1) 

^/i ^h-l 

i^,,H = -M-^ W^S,,,(0) := Mu,,n+i + P^h+i , (3-2) 
eo,h = — Ao,h := eo,/i+i + Po h+i : (3-3) 

^h-l 

ei,h = ^ := ei,h+i + l^i^h+i , (3-4) 

^h-l Vh-1 

and e2,h — ^i,h- The beta functions appearing in the r.h.s. of flow equations are related, 
see (2.7), to the kernels Wm,n4>,iJ.i so that they are expressed by series in the running cou- 
pling constants admitting the bound (2.27). For the explicit expressions of the one-loop 
contributions to the beta function, see below. 



B. The flow of 



IJ.,h 



Let us assume that e = maxfc<o{|e^,fc|} is small, that Z^/Zh-i < e and C < f h < 1 

for a suitable constant C, for all h < 0. Under these assumptions, the flow of v^^h can be 
controlled by suitably choosing the counterterms f^; in fact, if i/^ is chosen as 



= - E ^'''^U . (3-5) 

fe=— oo 

then the effective coupling i/^^/j is 

h 

= - E , (3.6) 



from which one finds that i'ij,^h can be expressed by a series in {e^,fe}fc<o, starting at second 
order and with coefficients bounded uniformly in h. At lowest order, if /i < and setting 



(see Appendix C): 



;r = -(M-l) 



^o,h'"h'^ l^h - arctan^;, 



-(M-i; 



dt (2x(0-x'(0) 



(3.7) 



27r2 



arctan,^/i — arctan ^/i 



dt {2x{t)-x'{t)) . 



(3.8) 
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By the above equations we see that lowest order contributions to i/^ are positive, that is 
i/fj, can be interpreted as hare photon masses. By using the short memory property and 
symmetry considerations, one can also show that ji^ f^ — j3\^^ is a sum of graphs whose 
contributions are of the order 0(1 — Vh) or O(eo,/i — e.i,h)- 



C. The flow of Zh and Vh 



In this section we show that, under proper assumptions on the flow of the effective charges, 
the effective Fermi velocity tend to a limit value fg// = V-oo and that both Ve// — and 
Zf^^ vanish as /i — )> — oo with an anomalous power law. 

Let us assume that the effective charges tend to a line of fixed points: 

e^,h = e^,_oo + 0{e\v^^ - vn)) + ©(e^M^'^) , (3.9) 

with < ^ < 1 and e^^_oo = e + 0{e^)\ this is a remarkable property that will be proven 
order by order in perturbation theory by using Wis, see the following section. Moreover, 
let be fixed as in the previous subsection (under the proper inductive assumptions on 
and Vk). 

We start by studying the flow of the Fermi velocity. At lowest order (see Appendix C) , 
its beta function reads: 

^^'^'^ = (3.10) 
logM 



47r2 



ep./i^fe^ arctan^fe / 2 ^o.ft^fe ^ \ ^fe - arctan 



Note that if eo,fe := ei^jt, then the r.h.s. of (3.10) is strictly positive for all > and it 

vanishes quadratically in at — 0. The higher order contributions to have similar 

properties. This can be proved as follows: we observe that the beta function is a function 
of the renormalized couplings and of the Fermi velocities on scales > h, i.e.: 

Ph = Ph{{ieo,k,ei,k,e2,k),{jyo,k,J^i,k,J^2,k),Vk}^^f) . (3.11) 

We can rewrite as /3l''^^ + /3l'^ + ^l'^ + /3l'^, with: 

^r^' = /3^({(eo,fc,eo,fc,eo,fc),(z/o,fe,z^o,fc,'^o,fc),l}fc>J , (3.12) 

l^h'^ = /3h({(eo,fc,eo,fc,eo,fc), (z/o,fe,t'o,fc,t'o,fe),^fe}fe>ft) ~ 

-Ph{{i^o,k, eo,fc, eo,A;), {i^o,k, ^o,k, ^o,k)i ^] k>h) ' 
Ph'^ ^ /3l{{{eo,k,eo,k,eo,k),{^o,k,i^i,k,'^2,k),Vk}^^,) - 

-f^h{{i^O,k, eQ,k, eo,k), {j^Q,k,J^Q,k,J^O,k),Vk}f,^f^) , 

Ph^ = ^h({(eo,fe,ei,ik,e2,ifc), (i^o,ik,i^i,fe,i^2,ifc),^^it}jt>;j) - 
-^h({(eo,ik,eo,ik,eo,fe), {i'o,k,vi,k,V2,k),Vk] . 
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By relativistic invariance it follows that ^'^'^^^ — and by the short memory property (see 
discussion after Eq.(2.32)) we get: 

= 0{el^{l - vn)) , /3f = 0{el,{uo,H - ^i,.)) , = O(eo,.(eo,, - e,,,)) . (3.13) 

Using (3.6) and an argument similar to the one leading to Eq.(3.13), we also find that 
^o,h ~ ^i,h can be written as a sum of contributions of order eo,/i(eo,fe — Ci^h) a-nd of order 
Bq ^(1 — Vh)- Therefore, we can write: 



Vh-i ^ ^ ^ ^ogM 
Vh 47r2 



^e'{l - Vh){l + A'J + le(l + S;)(eo,, - e^,,) 



(3.14) 



where: the numerical coefficients are obtained from the explicit lowest order computation 
(3.10); A'f^ is a sum of contributions that are finite at all orders in the effective couplings, 
which are either of order two or more in the effective charges, or vanishing at = 1; 
similarly, 5^ is a sum of contributions that arc finite at all orders in the effective couplings, 
which are of order two or more in the effective charges. From (3.14) it is apparent that Vh 
tends as /i ^ — oo to a limit value 

Veff = 1 + ^(eo,-oo - ei,_oo)(l + C'_J (3.15) 

with C'_ao ^ sum of contributions that arc finite at all orders in the effective couplings, which 
arc of order two or more in the effective charges. The fixed point (3.15) is found simply by 
requiring that in the limit h — oo the argument of the square brakets in (3.14) vanishes. 

Using Eq.(3.9), we find that the expression in square brackets in the r.h.s. of (3.14) can 
be rewritten as (8e^/5)(ve// — Vh + -R'/j)(l + A'l), where: (i) A'^ is a sum of contributions 
that are finite at all orders in the effective couplings, which are either of order two or more 
in the effective charges, or vanishing at — v^ff] (ii) R'^, is a sum of contributions that are 
finite at all orders in the effective couplings, which are of order two or more in the effective 
charges and are bounded at all orders by M^^, for some < ^ < 1. Therefore, (3.14) can 
be rewritten as 

Veff - VU-, = {veff - Vh) (l - ^ /^'^ J^!; ^^^^ log Mg(l + Ai)) , (3.16) 

from which, using the fact that = 0(e^M^'*), we get that there exist two positive con- 
stants Ci,C2 such that [35]: 

^^j^hf, < Veff - Vh ^ ^^^hf, ^ 
Veff - V 

with ^ 

r? = -logM [l-^;e//logM^(l + A':j] ; (3.18) 

at lowest order, Eq.(3.18) gives fj^'^^ = 2e^/(57r^). 
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Similarly CiM'^^ < < C2M^^ for two suitable positive constants Ci,C2, with rj — 
lim/i_;._oo logjv^(l + zo,h)', at lowest order we find (see Appendix C): 



,(2) _ logM 2 



2 „.-2\^h 



arctan 



(3.19) 



(2) 



127r^ 



so that 1] 

Before we conclude this Section, let us briefiy comment about the relation between Z^, 
Vh and the functions ^(k) and f (k) appearing in the main result, see (1.6). If |k| = M^, 
we define Z{k) = and f (k) = v^, for general |k| < 1, we let Z(h.) and f (k) be smooth 
interpolations of these sequences. Of course, we can choose these interpolations in such a 
way that, if M'* < |k| < M''+\ 





< 


Zh+1 ^ 






Zh 



v{k) - Vh 



Veff - Vh 



< 



Vh+1 - Vh 



Veff - Vh 



0{r)\ogM) , 
= 0(?7logM) 



(3.20) 



therefore, we can replace in the leading part of the 2-point Schwinger function (2.33) the 
wave function renormalization Zj and the effective Fermi velocity Vj by Z{k.) and f (k), 
provided that the correction term B{k.) in (2.33) is replaced by a quantity -B(k) defined so 
to take into account higher order corrections satisfying the bounds (3.20). This leads to the 
main result Eq.(1.6). 



D. The flow of the eff^ective charges at lowest order 

The physical behavior of the system is driven by the fiow of Cf^^h] hi the following section, 
by using a WI relating the three- and two-point functions, we will show that e^^h remain 
close to their initial values for all scales h < 1 and lim/i_^_oo e^^/j = e^^_oo = e + F^, where 

can be expressed as series in the renormalized couplings starting at third order in the 
effective charges. In perturbation theory, this fact follows from non-trivial cancellations 
that are present at all orders. For illustrative purposes, here we perform the lowest order 
computation in non-renormalized perturbation theory, in the presence of an infrared cutoff 
on the bosonic propagator xih,o]ip) — Xo(p) ~Xo(-^~'*p); at this lowest order, such a "naive" 
computation gives the same result as the renormalized one; for the full computation, see 
next section and Appendix D. If 

(7o, 7i, 72) := (70, ^^7l, ^^72) , (^0, h, h) := {ko, vki, w/ca) , (3.21) 
the effective charges on scale h at third order are given by: 

(eS - e)7. = ^e' J ^ [iJ^'^ {k)zj,9^^'^ {k)j^ + j.d.J^'^ (k)^.] (3.22) 
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where the first term in square brackets is the vertex renormahzation, while the second term 
is due to the wave function and velocity renormalizations. Note that both integrals are well 
defined in the ultraviolet (thanks to the presence of an ultraviolet cutofi^ in the propagators), 
while for h — )> — oo they are logarithmically divergent in the infrared. However, a remarkable 
cancellation takes place between the two integrals; in fact: 

^(^°)(k)z7,^(^°)(k) + a,^^(^°nk)= (3.23) 

9k Xo(k) .11 
= "-^ + Xo(k) (xo(k) - 1) ^iT;.^ , 

with Y '■— kojo + vk ■ J, so that 

C = e + ie' l ^l47„^xi(k)XM|(k) + 0(e'(M-l)) (3.24) 

Notice that the cancellation does not depend on the presence of the bosonic IR cutoff; this 
fact will play an important role in the analysis at all orders of the fiow of the effective 
charges, see next Section. An explicit computation of (3.24) says that, at third order in e, 

egL = e + eaf , (3.25) 

where 

= ^(2 - ( ^°""f"^° ) + 0(e\M - 1)) , (3.26) 

' - 4^^VA^ - ' 'i"""° ) + 0(AM - D) ; (3.27) 

the correction terms 0(e^(M — 1)) can be made as small as desired, by choosing < M — 1 <C 
1. Note that the two effective charges are different: 

e^loo - = -^^(^) + 0{e\M - 1)) , (3.28) 

where F{v) the function defined in (1.9). Combining (3.28) with (3.15) gives the last equation 
of (1.8) 

Of course, the one-loop computation that we just described docs not say much: if 
we could not guarantee that a similar cancellation takes place at all orders there would 
always be the possibility that higher orders produce a completely different behavior, e.g. a 
vanishing or diverging flow for e^^hi corresponding to completely different physical properties 
of the system. In order to obtain a control at all orders on e^^/j one needs to combine 
the multiscale evaluation of the effective potentials with Ward Identities. This is not a 
trivial task: Wilsonian RG methods are based on a multiscale momentum decomposition 
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which breaks the local gauge invariance, which Ward identities are based on. In Sec- 
tion 4 below, following a strategy recently proposed and developed in [6], we will prove (3.25). 

Remark. Note the unusual feature that e^^h ^ ei^^, an effect due to the presence of the 
momentum cut-off and the fact that v ^ 1. The discussion of this and previous sections 
can be repeated in the case that the bare interaction involves two different charges, cq and 
ei, describing the couplings of the photon field with the temporal and spatial components 
of the current. If e = (cq + ei)/2 and Cq — ei = O(e^), the conclusion is that fe// = 
1 — (e^/67r^)F(t;) + (5/6)(eo — ei)/e -|- O(e^) and it is of course possible to fine tune the 
bare parameters eo and ei in such a way that eo-oo = ei,_oo and Veff ~ 1. Note that, in 
a more realistic model for graphene, describing tight binding electrons on the honeycomb 
lattice coupled with a 3D photon field via a lattice gauge invariant coupling, one expects 
that eo-oo = ei _oo and v^ff = 1. 

4. WARD IDENTITIES 

In this section we prove that order by order in perturbation theory the effective charges 
e^^h remain close to their original values e^_o — e; moreover, we prove that asymptotically 
as /i — )■ — oo, eo,/i 7^ Ci^/j, see (3.25)-(3.28). The proof is based on a suitable combination 
of the RG methods described in the previous sections together with Ward Identities; even 
though the momentum regularization breaks the local gauge invariance needed to formally 
derive the Wis, we will be able, following the strategy of [6], to rigorously take into account 
the effects of cutoffs, and to control the corrections generated by their presence. 

As anticipated at the end of Section 2, we consider a sequence of models, to be called 
reference models in what follows, with different infrared bosonic cutoffs on scale h, i. e. with 
bosonic propagator given by: 

^[M](p) ^ Z^MiP) ^ XM](p)^Xo(p)-Xo(M-'^p) (4.1) 

(the idea of introducing an infrared cutoff only in the bosonic sector is borrowed from 
Adler and Bardecn [I], who used a similar regularization scheme in order to understand 
anomalies in quantum field theory). The generating functional yV[hfi] (^i 0) of the correlations 
of the reference model can be evaluated following an iterative procedure similar to the one 
described in Section 2 (see Appendix B for details), with the important difference that 
after the integration of the scale h we are left with a purely fermionic theory, which is 
superrenormalizable: in fact, setting m = in the formula for the scaling dimension (see 
lines following (2.27) and recall that for scales smaller than h the reference model has no 
bosonic lines) we recognize that the scaling dimension of this fermionic theory is 3 — 2n, 
which is always negative once that the two-legged subdiagrams have been renormalized, see 
[16]. Let us denote by {e^a\}k>h the effective couplings of the reference model; of course, if 
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k>h 

= e,,, , (4.2) 

where {e^^fe}fc<o are the running couphng constants of the original model. On the other hand, 
as proven in Appendix B, the vertex functions (j/i,p; V'k+pV'k)^ reference model with 

bosonic cutoff on scale h computed at external momenta k, k+p such that |k+p|, |k| ~ 

and IpI ^ M'^ are proportional to the charges e|^'^ = e^^/j, see (2.34); therefore, if we 
get informations on the vertex functions of the reference models, we automatically infer 
informations on the effective couplings of the original model. 

Such informations are provided by Ward Identities; by performing the change of variables 
ipx e*"'''^x, 'ipx — ^ e~'^°''^il)y. in the generating functional W[/i.o](J, 0) of the reference model 
and using that the Jacobian of this transformation is equal to 1, see [6], we get: 



J p(dV')P[^,o](dA)e-^'^'^^-(^"'"^^^'"^-^)'^-+^(^''^)+^(-^''^^"'") , (4.3) 



where P[hfl]{dA) is the gaussian integration with propagator (4.1) and, iiY — 7o^o + v^-k, 
the pseudo-differential operator D is defined by: 

f dk. e^^^ 
4(k)>o {^^r Xo(k) 

If we derive (4.3) with respect to a, (f) and (f) and then set a = = J = 0, we get the 
following identity: 



where 



and 



PmO'/^,p; V'k+pV'k)fe = (V'kV'k)/, - (V'k+pV'k+p)ft + ^h{K p) (4.4) 

/dk.' — — 
(V'k'+pC(k', p)V'k'; V'k+pV'k)^ (4.5) 



C(k, p) = ^K(xo(k)-^ _ i) _ ,(|^ + j^) (^;^,(k + p)-i - i) . (4.6) 

The correction term A/i(k, p) in (4.4) is due to the presence of the ultraviolet momentum 
cut-off, and it can be computed by following a strategy analogous to the one used to prove 
the vanishing of the beta function in one-dimensional Fermi systems [6] . We can write 

A/,(k, p) = a^p^,{j^,,p] V'k+pV'k)ft + -7^Ri.,h{^, P) , (4.7) 

where the correction i?^ /i(k, p) is dimensionally negligible with respect to the first term, 
see Appendix D. More precisely, in Appendix D it is shown that: (i) -R^^/j(k, p) can be 
written as a sum over trees with N endpoints of contributions i?^^''(k, p); (ii) it is possible 
to choose in such a way that, under the same conditions of Theorem 2.1 and if |k| = M^, 
|k + p| < and |p| < M^, 

|<,)(k,p)| < (const.)^(^)!M-^'^Mt£t . (4.8) 
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An explicit computation, see Appendix D, shows that at lowest order a^^^ is given by (3.26)- 
(3.27). 

Let us now show how to use the previous relations in order to derive bounds on the 
effective charges. Let us pick |k| = and |p| ^ M'*; by using Eqs.(2.33)-(2.34) and the 
fact that 

^W(k) _^W(k + p) = ^W(k + p)(^po7o + ^^/.-lP•7)5^'Hk) +P,r,{\^,^) , (4.9) 
with f^(k,p) = 0{\p\M-^^), we find that 

(^k?k)/^ - (V-k+pV^k+p)^ = i^g^^K"^ + p)(^Po7o + tVh-iP- Dg'^^K'^) + 

■^/t— 1 

+ ^(rVk,p)+f^(k,p)) , (4.10) 

Pn{jn,p\ V'k+pV'k)ft = -^—9^''\^ + P)(^eo,fePo7o + ivh-iei,hP- l)g^^\^) + 

+^r,(k,p), (4.11) 

with |r^(k, p)|, |f^(k, p)| expressed by sums over trees of order A?" > 2 of contributions 
r^^'*(k, p), f^^'*(k, p) bounded by (see Appendix B, formulas (B.13) and (B.16)) 

\rf \Kp)\ + |rf Hk,p)| < (const.)^£t(f j'M-^'^ . (4.12) 

Now, if we plug (4.7) into the Ward identity (4.4), and we use the relations (4.10)-(4.11), 
we get an identity that, computed at k = ko := (M'*, 0) and p = po := (p, 0), after taking 
the limit p — >■ 0, reduces to: 

^(1 - ao) = 1 + iM^"" [ro(ko, 0) + i?o,/^(ko, 0) - (1 - ao)ro(ko, 0)]7o 

= 1 + Ao,h , (4.13) 

with A^^h a sum of contributions associated to trees of order N > 2 bounded at the A^- 
th order by (const. )^(£/i)^(A/2)!, as it follows from the estimates on RQ,h-i^O)fo (note the 
crucial point that such estimate is proportional to {sh)^ rather than to (£_oo)^; this is the 
main reason why we chose to introduce the infrared cutoff on the bosonic propagator, see 
the end of Section 2E and the beginning of this section). Eq.(4.13) combined with (3.26) 
implies, as desired, that the effective charge eo./j remains close to eo,o = e at all orders in 
renormalized perturbation theory. Moreover, proceeding as in the derivation of Eq.(3.13), we 
find that |Ao,/i — ,-oo| = 0{e^{veff — Vh)) + 0{e^M^^), for some < < 1, from which we 
get Eq.(3.9) for /i = 0. Similarly, if ki := (0, M^, 0), we get: 

^(1-q;i) = h^M Vh-i- 

e Vh 

■ [fi(ki, 0) + i?i,,,(ki, 0) - (1 - ai)ri(ki, 0)]7i 
= 1 + A,^, (4.14) 
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with Ai^h a sum of contributions associated to trees of order N > 2 bounded at the A^-th 
order by (const. )^(£/i)^(A^/2)!, which imphcs that the ejfective charge ei^h remains close to 
ei^o = e at all orders in renormalized perturbation theory. Moreover, as in the = case, 
\M,h — ^1 -ool = 0{e'^{veff — Vh)) + 0{e^M^^), for some < ^ < 1, from which we get 
Eq.(3.9) for /i = 1. 

Equations (4.13) and (4.14) not only imply the boundedncss of the effective charges e^^^, 
but they also allow us to compute the difference eo,/j — ci^u, asymptotically as /i — )■ — oo, at 
all orders in renormalized perturbation theory. At lowest order, Cq^^ — e[^^ = e(Q;Q^'' — af''), 
as anticipated in previous section. 

5. CONCLUSIONS 

We considered an effective continuum model for the low energy physics of single-layer 
graphene, first introduced by Gonzalez et al. in [17]. We analyzed it by constructive Renor- 
malization Group methods, which have already been proved effective in the non perturha- 
tive study of several low-dimensional fermionic models, such as one-dimensional interacting 
fermions [6], or the Hubbard model on the honeycomb lattice [16]. While in the present 
case we are not able yet to prove the convergence of the renormalized expansion, we can 
prove that it is order by order finite, see Theorem 2.1 above. Note that, on the contrary, 
the power series expansion in the bare couplings is plagued by logarithmic divergences and, 
therefore, informations obtained from it by lowest order truncation are quite unreliable. In 
perspective, the proof of convergence of the renormalized expansion appears to be much 
more difficult than the one in [6] or [16], due to the simultaneous presence of bosons and 
fermions, but it should be feasible (by using determinant bounds for the fermionic sector 
and cluster expansion techniques for the boson sector). 

A key point of our analysis is the control at all orders of the flow of the effective couplings: 
this is obtained via Ward Identities relating three- and two-point functions, by using a 
technique developed in [6] for the analysis of Luttinger liquids, in cases where bosonization 
cannot be applied (like in the presence of an underlying lattice or of non-linear bands) . The 
Ward Identities have corrections with respect to the formal ones, due to the presence of a 
fermionic ultraviolet cut-off. Remarkably, these corrections can be rigorously bounded at all 
orders in renormalized perturbation theory (see Section 4). 

Several questions remain to be understood. First of all, the effective model we considered 
is clearly not fundamental: a more realistic model for graphene should be obtained by 
considering electrons on the honeycomb lattice coupled to an electromagnetic field living in 
the 3D continuum. We believe that a Renormalization Group analysis, similar to the one 
we performed here, is possible also for the lattice model, by combining the techniques and 
results of [16] with those of the present paper; we expect that the lattice model is asymptotic 
to the continuum one considered here, provided that the bare parameters of the continuum 
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model are properly tuned. Another important open problem is to understand the behavior 
of the system in the case of static Coulomb interactions; this case can be obtained by taking 
the limit c — )■ oo together with a proper rescaling of the electronic charge in the model with 
retarded interactions. However, as discussed in the Remark at the end of Section 2D, the 
static case seems to be much more subtle than the one considered in this paper, since it 
apparently requires cancellations even to prove renormalizability of the theory at all orders. 
We plan to come back to this case in a future publication. 
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Appendix A: Symmetries 

In this Appendix we prove formulas (2.6) and (2.7); to do this, we exploit suitable sym- 
metry transformations. We use the following explicit representation of the euclidean gamma 
matrices: 

where ui = and (72 = . It is also useful to dehne: 73 = I , with 

yi oy \^ ^ / \-«f^3 J 

(Js = , and the corresponding fifth gamma matrix 

75 = 7o7i7273 = , (A.2) 

which anticommutes with all the other gamma matrices: {7^, 75} = 0, V/i = 0, . . . , 3. Finally, 
given u G {+, — }, we define the chiral projector := (1 + u'~f^)/2. 

It is straightforward to check that both the gaussian integrations P{dilj), P{dA) and the 
interaction V{A,iIj) are invariant under the following symmetry transformations, which are 
preserved by the multiscale integration: 

(1) Chirahty: P^jipk P^^ipk, ip^P-uj V^k-P-^e"^*"'^ , with e M independent of k. 



(2) Spatial rotations: — >■ e4['>'i''>'2]'0^[i,2]j^, ■0k — >■ ijj ^[i,2]^^e ll'^^'^^^ and A/^^p — >■ 



[4'''^^..#,^iJ,,with 



/l \ 

cos^ — sin^ 
^0 sin 6* cos 6* J 



(A.3) 
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The invariance of the model under (2) is a simple consequence of the fact that 

g-f [71,72] (-^^^^^^^^j gf [71,72] ^ (70, -)/^ cos 61 - 72 sin 61, 72 cos 61 + 7i sin 61) . (A.4) 

(3) Complex conjugation: ip^ {—il2)4'-'k, "^k ~^ V^-k(^72)i ^^.k ~^ — ^Ai,-k and k k,*, 
where k is a generic constant appearing in P{dip), P{dA) and/or in V{A,ip). 

(4.a) Horizontal reflections: ifj]^ («737i)'0k, '^k ^ and >l^,p (-l)''A^,p, 

where k = {ko, —ki, /C2). 

(4.b) Vertical reflections: ■^k — >■ {—il2)'4>i^: V'k ^ i^xS^l'i) and ^^,p — >■ (— l)'''''2 74ju^p, where 
k = (fco, fci, -A;2). 

(5) Particle - hole: V'k ^ (-7072)^^1^ V'k ^ "0^7270, and .4^,p ^ (-l)^-^^'M^,p, where 
k = (/co, -A;). 

(6) Inversion: Vk 7o73V'k> V'k ^ V'k737o and A^,k (-l)'^'''M^,p, where k = {-ko, k). 

In addition to the previous symmetries, if v = c = 1 the theory has an additional space-time 
invariance, namely: 

(7) Relativistic invariance: — >■ ^^^'^^'^^'^i^ j^yo^]]^-, V'k ~^ V'^io-iik^"^'''^''''''^^ and A^^-p — >■ 
[^eA.^^-ip]^, with 

/ cos ^ — sin ^ 0^ 



p[o,i] _ 



sine cose I . (A.5) 
V 1^ 



The invariance of the model under (7) is a simple consequence of the remark that 

g-f [70,71] (^^^ ry^^ 1'^°''^^] = (70 cos ^ — 7i sin e, 71 cos ^ + 70 sin e, 72) . (A.6) 

It is now straightforward to check that these symmetries imply (2.6), (2.7). In fact, the first 
two identities in the first line of (2.6) and the second identity in the second line of (2.6) 
easily follow from (4.a)-|-(4.b) + (6). Using (4.a)-|-(4.b)-|-(6) we also find that 

Mt,M = , (A.7) 

while, from (2)-|-(3), we get 

W^iSi,i(0) = ^i,o,2,2(0) , W^i,o!,,,(0) e R , (A.8) 

which imply the first identity in the second line of (2.6) (notice that, if v = c = 1, from (7) 
we also get that wi%^,{0) = <'o!i,i(0))- 
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Let us now consider the combination ■0k^o,?(k)'0k- Using the fact that 
{/, 75, {7j}o<j<3, {7j75}o<j<3, {7ji7j2}o<ji<j2<3} is a complete basis for the space of complex 
4x4 matrices, we can rewrite it as: 

3 

V^k<?(k)V'k = ^^{co{^)I+Yl H(^hM,(^hM+ E 4^'^(k)7..7,.+C5(k)75}V'k ■ 

j=0 0<ji<j2<5 

(A.9) 

Now, using the invariance under (1), we find that, e.g., 

V^kCo(k)/V'k = V^k^a.Co(k)7P,V'k = e-2-<^ ^ V^k^,Co(k)7P,V'k , (A.IO) 

for all e M, which implies that Co(k) = 0; similarly, using the invariance under (1) and 
the fact that [75,^0;] — [ijilji^PiA = 0, VO < ji < j2 < 3, we find that C5(k) = and 
= 0, VO < ji < j2 < 3. Therefore, 

3 

/ (^^•'<i^W^>< = / E(^V^k[ci(0)7, + 4(0)7,-75]V'k, (A.ll) 

3 

/ (^^•^k9,<';)(0)^k = pj ^?k[k9kci(0)7, + kak4(0)7,75]^k . (A.12) 

Let us first look at (A.ll). Using the invariance under (4. a), we find that [c{(0)7j + 
4(0)7^75] = 737i[ci(0)7i + fi5(0)7i75]7i73, which implies that c[{0) = c?(0) = 0^5(0) = 
^15(0) = 0. Using (2), we find that also cl{0) = cl^{0) = 0; finally, using (6), we find that 
c5(0) = 4(0) = 0. This concludes the proof of the third identity in the first fine of (2.6). 

Let us now look at (A.12). The terms proprtional to ko in the r.h.s. of (A.12) are invariant 
under (2)+(4.a)+(4.b), which implies that dkAi^) = dkA{0) = dk^cHO) = 9^^4(0) = 0. 
The terms proportional to ki are invariant (4.b)+(6), while the terms proportional to k2 
are invariant (4.a)+(6); combining these transformations with (2), we find that 9fe^c°(0) = 

dkAiO) = dkAiO) = 9fei4(0) = 0, that 9feC?(0) = dkA{0) = d^AiO) = ^^2^^15(0) = 0, 
and that dk^cKO) — dk^cKO). Therefore, 

^V^kkak<? (O)V^k = j ^V^k [«o^o7o + ai^ • 7] V'k , (A.13) 

for two suitable constants ao,ai. Using the invariance under (3), we find that ag = izQ^h 
and ai = izi^h, with Zfj_^h G M, which concludes the proof of the first line of (2.7) (of course, 
if V = c = 1, then from (7) we also get that zo,h = that is the speed of light is not 
renormalized) . 

A completely analogous discussion can be repeated for the second line of (2.7), but we 
will not belabor the details here. 
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Appendix B: Multiscale integration for the correlation functions 

The multiscale integration used to compute the partition function W(0,0), described in 
Section 2, can be suitably modified in order to compute the two and three-point correlation 
functions in the reference model with bosonic infrared cutoff on scale h, see (4.1). We start 
by rewriting the two and three point Schwinger functions in the following way: 

where W[/i*,o](>^, 0) is the generating function of the reference model. The two point 
Schwinger function (V'kV'k) appearing in our main result is obtained as lim/j«_^_oo (V'kV'k)/!*- 
In order to compute W[h*,Q\{J-i(t>)-i we proceed in a way analogous to the one described 
in Section 2. We iteratively integrate the fields -0*^°^ ^4*^°', . . ., ■0^'^+^), ^l^'^"'"^^ . . ., and after 
the integration of the first \h\ scales we are left with a functional integral similar to (2.2), 
but now involving new terms depending on J, 0. Let us first consider the case h > h*; the 
regime h < h* will be discussed later. 

Case h > h*. We want to inductively prove that 



e 



^ih*,oM<l>) = (B.2) 



,|A|£;;,+5(>'^)(J,^ 



P((iV'(^''y[,.,o](fi^^-'^)e^*''^^*'''+''"'''^^*'"^- 



+ 



where: S''-^'' {J, (f)) is independent of {A,iIj), W^^ contains terms explicitly depending on 
{A, -0) and of order > 2 in 0, while is given by: 

4'^)(AV%V',0) = I ^[0k[g('^+^)(k)]Vk + V^k^3^'+'Hk)0k 

+^V^'^\A,^,i;)Gf^'\k)<j>^] . (B.3) 
Moreover, the functions Ga, Q^^\ G^^^ are defined by the following relations: 

i=h 

gW(k) := g('^+i)(k) - iZ^z.^^k.^.G^^-^'Xk) , (B.4) 
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with g(^)(k) = 1, ^(^^(k) = 0. Note that, if k is in the support of g^''\\s.), 



(B.5) 



that is \\Q^''\'k) - 1|| < (const.) 4 and ||Gj!'^(k)|| < (const.) Z^^M"'*. Moreover, by the 
compact support properties of w^^*'^\p), GA,At(p) = for all |p| < M^* . 

In order to prove (B.2)-(B.4) by induction, let us first check them at the first step. The 
generating functional of the correlations is defined as (see (1. !)-(!. 2)) 



(B.6) 



P(dV^(^°))P[,.,oi(^^^-°^)e' 



.(<()) 



(<o) a(<0) 



which, under the change of variables 

can be rewritten in the form (B.2), with 

= , e^.S^^'^) = iy,{J„ J,) + uliJ,, J^) , 

WP^O, V(°) = F, M°)(^ + G^J,^,0) = 5(0,0) 



(B.7) 



(B.8) 



Let us now assume that (B.2)-(B.4) are valid at scales > h, and let us prove that the 
inductive assumption is reproduced at scale h — 1. We proceed as in Section 2; first, we 
renormalize the free measure by reabsorbing into P(di/j^-^^) the term exp{C,i,V^'^^}, see (2.8)- 
(2.11), and then we rescale the fields as in (2.12). Similarly, in the definition of -B^'*'', 
Eq.(B.3), we rewrite V*^''^ = C^V^^^ + V*-'*-', combine the terms proportional to >C^V^'*-' with 
those proportional to Q^'^^^\ and rewrite 



dk r_ _ 

^>k[g(''Hk)]Vk + ^kQ^'^(k)0k 



+ 



(27r)3 
rik 



(27r)3 
d 



0k[Gr'^(k)]^-^v('^)(A, yz;::^^) + 



d 



with Q^'^^ defined by (B.4). Finally, we rescale by defining 

wi^\A + GaJ, ^fZ^xi^) := wf{A + G^J, ^^i^) , 
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and perform the integration on scale h: 

_^B('')(A(<'')+GAJ,v/^V'(^''))+l^i''' = 
^ g|A|fih+5(''-i)(J,</.)+V(''-i)(A(<''-i)+GAJ,A/2;rTV'(^''-i') . 

_^B(''-i)(A(<''-i)+G^J,^/i^V(^''-i))+<''-'^ ^ 

where S^^^^\J, (p) contains terms depending on (J, 0) but independent of ^4^-'*"^^ 
Defining S^-^'^"* := S'^^'^^ + S^-'^\ we immediately see that the inductive assumption is 
reproduced on scale h — 1. 

Case h < h*. For scales smaller than h*, there are no more bosonic fields to be integrated 
out, and we are left with a purely fermionic theory, with scaling dimensions 3 — 2n, 2n being 
the number of external fermionic legs, see Theorem 2.1 and following lines. Therefore, once 
that the two-legged subdiagrams have been renormalized and step by step reabsorbed into 
the free fermionic measure, we are left with a superrenormalizable theory, as in [16]. In par- 
ticular, the four fermions interaction is irrelevant, while the wave function renormalization 
and the Fermi velocity are modified by a finite amount with respect to their values at h*; 
that is, if Eh* = maxfe>/,.{|e^,fc|, Wn,k\}' 

Zh^Zh*{l + 0{el)) , Vh^Vh^{l + 0{el)) . (B.9) 

Tree expansion for the 2-point function. As for the partition function (see section 2), the 
kernels of the effective potentials produced by the multiscale integration of yV[h*fi]{J, 0) can 
be represented as sums over trees, which in turn can be evaluated as sums over Feynman 
graphs. Let us consider first the expansion for the 2-point Schwinger function. After having 
taken functional derivatives with respect to 0k, 0k ^^^d after having set J = = 0, we get 
an expansion in terms of a new class of trees r G T^'^ ^, with k e {—oo, —1] the scale of the 
root and h > k; these trees are similar to the ones described in section 2, up to the following 
differences. 

(1) There are N + 2 end-points and two of them, called vi,V2, are special and, respectively, 
correspond to [g('^''i-i)(k)] V?'''""'^ or to V^if'''''"'^Q('*''2-i)(k). 

(2) The first vertex whose cluster contains both vi, V2, denoted by v, is on scale h. No TZ 
operation is associated to the vertices on the line joining v to the root. 

(3) There are no lines external to the cluster corresponding to the root. 

(4) There are no bosonic lines external to clusters on scale h < h*. 

In terms of the new trees, we can expand the 2-point Schwinger function as: 

(^kV^k).. = E[^?(k)ry^^^''nk) + E E E E (b.io) 
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where /ik < is the integer such that M^^ < |k| < M'*''+^, and <S2(r;k) is defined in a 
way similar to V^'^^t) in (2.23), modulo the modifications described in items (l)-(4) above. 
Using the bounds described immediately after (B.5), which arc valid for k belonging to the 
support of g^'^^'k), and proceeding as in Section 2D, we get bounds on (S2(T;k), which are 
the analogues of Theorem 2.1: 

E E E l|52(r;k)||<(const.r.t*(^)!^. (B.ll) 

h=—oo k=—oo T-cTlCi*) 

^'k,h,N 

Notice that the result (B.IO) and the bound (B.ll) are true for any k such that |k| > M'^* ; 
being the bound (B.ll) uniform in h*, our result on the two point function (1.6) and (2.33) 
is obtained by fixing k and taking the hmit h* — > — oo in (B.IO). 

In order to understand (B.ll), it is enough to notice that, as far as dimensional bounds 
are concerned, the vertices Vi and V2 play the role of two u vertices with an external line (the 
line) and an extra Z^^^^M~^^ factor each. Moreover, since the vertices on the path Vr,v 
connecting the root with v are not associated with any TZ operation, we need to multiply the 
value of the tree r e 7^^^*J^ by M(^/2){;i-fc)^(i/2)(fe-fc)^ ^^^^ exploit the factor M(V2)(fe-/i) 

in order to renormalize all the clusters in Vr,v- Therefore, 



E E E ll'52(r;k)||<(const.)^(-)! 



h=—oo k=—oo T-cTl('i*) 

^^k^h-h^j^{l/2){h-k)j^-2hy, (B.12) 
h<hk k<h 

where: the factor is due to the fact that graphs associated to the trees r e T^'^ ]^ 
have two external lines; the factor M^~^^ is given by the product of the two short memory 
factors associated to the two paths connecting v with vi and V2, respectively; the "bad" 
factor M(i/2)(?^-fc) is the price to pay to renormalize the vertices in Vr,v] the ^ and the 
last M"^^'' are due to the fact that Vi, V2 behave dimensionally as u vertices times an extra 
Zf^^ factor. Performing the summation over k and h in (B.12), we get (B.ll). Note 

also that, if k is on scale h^. — h*, then the derivatives of ||52(r;k)|| can be dimensionally 
bounded as 

E E E P."'52(r;k)||<(const.)^4(y)!^%;;;^, (B.13) 

h=—oo k=—oo tcT^^*^ 
^ 'k,h,N 

from which the bound on f^^^(k, p) stated in (4.12) immediately follows. 

Tree expansion for the ?> -point function. Let us pick |k| = M'^*, |k + p| < M^* and 
|p| <^ M^* , which is the condition that we need in order to apply Ward Identities in the 
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form described in Section 4. In this case, the expansion of 3-point function (jju,p; ■0k+p'0k)/i* 
is very similar to the one just described for the 2-point function. The result can be written 
in the form 

0,,p; ^k+p?k),. = [af "'^k + p)]t7,^?(''*) (k)Q('^') (k) + 

+ E E E '53(r;k,p) (B.14) 

N>1, k<h, t-gt'''*' 

h<h* hv^>h* k,h,hv3,N 

where Tr'l I is a new class of trees, with k < the scale of the root, similar to the 

k,h,h^^,N ' ' 

trees in 7^*^^ up to the fact that they have N -\-?> endpoints rather than N + 2 (see item 
(1) in the list preceding (B.IO)); three of them are special: Vi and t'2 arc associated to 
the same contributions described in item (1) above, while is associated to a contribution 
Zh^^-i (e^,/tii3/e)j^"p^"^''-M^^3(z/^ ;^_ye)A^^p, with the vertex immediately preceding on 
r (which the endpoint fa is attached to) and h^,^ > h*. The value of the tree, S^It; k, p), is 
defined in a way similar to S2{t; k), modulo the modifications described above. S^It; k, p) 
admits bounds analogous to (B.11)-(B.12); recalling that |k| = M^* , |k + p| < M^* and 
IpI < M''*, we find: 

E E E E l|53(r;k,p)||<(const.)-(f)!4^- 

h=-oo k=-oo hv3=h*+l ^^-rih*) 

. ^ jv^(V2)(fe-/i)7^^ft-/»*7v^(i/2)(/»*-fe„3)^-2fe* ^ (B.15) 

h<h* 
k<h 
hv3>h* 

where: M^^/^^^*^"'*) is the short memory factor associated to the path between the root and 
v; M^"^* is the product of the two short memory factors associated to the paths connecting 
V with vi and V2, respectively; M*^^/^)('**~'*''3) is the short memory factor associated to a path 
between h* and v^; M'"^^* /Zh*_i is the product of two factors M~'^^Z^^[\ associated to the 
vertices Vi and V2 (see the discussion following (B.ll) and recall that in this case h]^ = h*). 
We remark that in this case, contrary to the case of the 2-point function, the fact that there 
is no TZ operator acting on the vertices on the path between the root and v does not create 
any problem, since those vertices arc automatically irrelevant (they behave as vertices with 
at least 5 external lines, i.e., J, (f), (f) and at least two fermionic lines) and, therefore, TZ = 1 
on them. Note also that the vertices of type J(j)ip^ which have an TZ operator acting on, can 
only be on scale /i* — 1 or h* (by conservation of momentum) and, therefore, the action of 
the TZ operator on such vertices automatically gives the usual dimensional gain of the form 
const. M^"'^^' . Performing the summations over k^h^hy^ in (B.15), we find the analogue of 
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(B.ll): 



E E E E ll<S3(r;k,p)||< (const.)- 



h=-<x> k=-oo K3=h*+1 i-gTiCi") 

k,h,hv^ ,N 



Zh*- 



(B.16) 



from which the bound on ri-^(k,p) stated in (4.12). 

Appendix C: Lowest order computations 

In this Appendix we reproduce the details of the second order computations leading to 
(3.7), (3.8), (3.19), (3.10). 



1. Computation of 

By definition, see (2.7) and (3.1), = Zq^h — — ^7o^ito^o,'i\o)- one-loop, defining 
eo,/i = eo,/i and ei^h = Vh-iei,h, we find: 



^,(2) _ ^(2) 



4S = -noet.« / (^»P.(^)7.9"""(P)7. + (C.l) 



■ -2 / Zfi+l 



dp a (fh+2{v 



J J (27r)3^^«V 2|p 
Zh+i [dp (fh+i{p 

Co 



noe^,,+2 y (27r)3''^»l 2|p 



Using inductively the beta function equations for Z/j+i, v^j+i, 6^,^+2, and neglecting higher 
order terms, we can rewrite (C.l) as 

Po ^7m7o7m 



(2) _ . -2 1 /■ Pi 

^o,/.-noe^,/.+i2 y (27r)3|p 



(A+1(P) - |p|/ft+l(p))(A+l(p) + A+2(P)) 

+(A+2(p) - |p|/;:+2(p))A+i(p)" 



(C.2) 



Passing to radial coordinates, p = p(cos ^, sin ^ cos 99, sin 6' sin*/)), and using the fact that 
/ dp{f'f^^Jh+i + fh+ifh+2 + f'h+2fh+i) = 0, we find: 



»1 

dcos9 



87r2 



dp 
P 



cos 



cos2 6' + w^sin 9- 



(C.3) 
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The integral over the radial coordinate p can be computed by using the definition (2.1): 

r-(A%i + 2W.+2)= (C.4) 

JO P 

= r -[2(x(p) - x(Mp)) - ix'ip) - x'iMp))] = hm ^ = logM . 

Jo P ^^^Je P 

Finally, an exphcit evaluation of the integral over dcos9 leads to (3.19). 



(2) 

2. Computation of zl ^ 

By definition, see formulas (2.7) and (3.1), ^1'^^^ — zf\ — VhZ^l, with zi^h — 
—'i'lidkiW^'^i (0) . At second order, proceeding as in the derivation of (C.2), we find: 



(2) _ . -2 ^ 

^i,h — '^7ie^,/i+i^ • 



dp pI ilt,Vh-ii-i^ 



{2i:fWpl + vl\p\^ 

+(A+2(p) - |p|/;:+2(p))A+i(p) 



— {fl+l + 2//i+i//j+2) 
P 



(A+i(p) - iPi/;:+i(p))(A+i(p) + A+2(p)) + 



An explicit evaluation of the integral leads to 



d cos 6 



sill" 



e 



-1 



cos^ 9 + vl sin^ 9 



.(2) 



^1 ^ — ^0,hV. 

which, combined with ~ "'^'^^ 



^ log M / arctan — arctan 



S leads to (3.19). 



(C.5) 



(C.6) 



3. Computation of d'^'^i^^ 

By definition, see (2.6) and (3.1), ^^^^^ = -M-^+^W^'j;^i^l{0) - Mu^^h- At second order, 
we find: 



'^''^ 2 J (27r)3 



Tr(7,5^'^Hp)7M^^'Hp)) 



+ 



(C.7) 



Using inductively the beta function equations for e^^/j, and neglecting higher 

orders, we can rewrite (C.7) as 

«-.(2) _ _^j^-h+12 f dp A(p)' + 2/fe(p)/fe+i(p) 2 2 2 

W2)_ 9A^-/»+lp2 /" C^p /fe(p)' + 2/fe(p)/fe+i(p) 2 .p.^ 
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where we used that Tr (7^70,7^7^^) — —A ii /i a and 4 otherwise; passing to radial coordi- 
nates we find 



Po,h — 



3^^,(2) 
'l,h 



-2 

/: 

-2 

/: 



M- 



-h+i 2 



dcosO 



- cos^ d + vl sin 9 
[cos^O + f^sin^6')2 



dp {fl + Vhh+i) 



dcosO 



cos^ 9 



(C.9) 



_i (cos2 6' + v^sin''6') 

The integral over the radial coordinate p can be rewritten as, using the definition (2.1): 

poo 

dp {fl + Vhh+i) = M^-\M - 1) / dp (2x(p) - x\p)) ■ (CIO) 
Finally, an explicit evaluation of the integral over dcosO leads to (3.7). 

Appendix D: Multiscale integration of the correction term to the WI 

In this Appendix we prove (4.7) and the bound (4.8). We assume that h = h*, |k| = M'* 
and |p| -C M^. We start by rewriting 



-^R^ h* (k, p) = — - — = 

Zh ' dJpd(j)^^pd(j)k 



W[Ml(-/,0)L. 



J=(A=0 



(D.l) 



with 'W[h,Q]{J , (p) defined as: 

ey^lHM(l<t>) J P{d'iP)P[h,o]{dA)~ 



,V{A,^)+B{J,4,) 



(D.2) 



and 



+ 



+ 



dk 



(D.3) 



The main difference with respect to the generating functional of the correlation functions 
is the presence of the correction term proportional to C(k, p), see (4.6) for a definition. 
Eq.(D.2) can again be studied by RG methods, see [6] for further details. A crucial role is 
played by the properties of the function C(k, p); it is easy to verify that 



^«(k + p)C(k,p)^(^)(k) 



(D.4) 
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is non vanishing only if at least one of the indices i,j is equal to 0; moreover, when it is 
nonvanishing, it is dimensionally bounded from above by (const.) |p|M~'~^. 
We start by integrating the scale 0, and we find: 

gW[;,,o](J,<^) ^ g|A|B_i+5(^-i)(J,<^) . 

where .S^-"^) collects the terms depending on J, but independent of A, t/j, and 

{A, ^) = {A, ij, <P) + B^-'^ {A, ^) + W^'^ , (D.6) 

with: B[j ^\A, ip, (p) linear in J and independent of 0; ^\A, tp) given by (B.3); wj^^^ the 
rest, which is at least quadratic in {J ■,(()). With respect to the computation of W[h,o](<^5 0), 
we now have new marginal terms of the form Jipip, which are contained in 0) 
and need to be renormalized. Let us simbolically represent by W^m,rt''({kj}, {qj}, p) the 
generic non-trivial kernel appearing in i5j ^''(A,-?/', 0); m is the number of bosonic external 
lines (of either J or A type) while 2n is the number of %p fields; {kj}, {qj} are respectively 
the fermionic/bosonic momenta and p is the momentum flowing through J. As usual, 
these new kernels can be represented as sums over Feynman graphs. The J external line 
can be attached to a simple vertex, corresponding to the monomial —ctfj.PiiJpjl^p\ or to a 
"thick" vertex, representing Jp'0k^pC(k, p)'0k (the "small circle" associated to the vertex 
represents the matrix kernel C(k, p), see Fig. 6). Let us denote by Wm,n''^ the contribution 
to Wm,n coming from graphs with the J line attached to a thick vertex, see Figure 6. By the 




V 



FIG. 6: Schematic representation of the expansion for Wq ^ ; the small circle represents C(k, p). 

properties of the C(k, p) function, see [6] for details, it follows that wiiTn^i^j}) {li}) P) — ■ 
PiiWm,n,fM{{ki}, {qi},p), with Wrn,n,^i dimcnsionally bounded as an Aipip kernel, uniformly in 
p. We define the action of the TZ = 1 — C operator on Wm,n,iJ. in a way similar to (2.5)-(2.7). 
In particular, /^W^Q^^^j (k, p) := Wq'^^^{0,0) and, by symmetry, 

Z_i JpV^k£<-;>k+p = -Z_2 Jpa^,-4^-^) , (D.7) 

for a real constant which is by definition the effective a-coupling on scale —1. Note 

that the last two graphs in Fig. 6 do not contribute to a^u^-i simply because they are one- 
particle reducible and, therefore, they are vanishing at zero external momenta. 
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We now iterate the same procedure, and step by step the local parts of the kernels of type 
are collected together to form a new running coupling constant, a^^^; in order to show 
that R^^h is dimensionally negligible as /i — )■ — oo, we need to show that it is possible to fix 
the initial data = o in such a way that a^^h goes exponentially to zero as 77 — >■ —00, 
which is proved in the following. 

The flow of a^^k ■ The new marginal running coupling constants evolve according to 
the flow equation: a^^k-i = ot^^k + f^'^^ki where a^^ = are the counterterms appearing in 
the bare interaction (D.3). The beta function /S^^^ can be split as 



3 a 



+ 



oa,2 



(D.8) 



where collects the contributions independent of a^,^' (which, therefore, are associated 
to graphs with the J external line emerging from the thick vertex representing C(k, p)), and 
P'^'l collects the terms from graphs with one vertex of type a^^k' ioi some k' > k. It is crucial 
to recall that by the properties of C(k, p), the graphs contributing to jS^'l have at least one 
propagator on scale or —1; by the short memory property, this means that they can be 
dimensionally bounded by (const. )e^M''^, for any < ^ < 1. Similarly, the contributions to 
P'^'l associated to graphs with at least one vertex of type a^^k' for some k' > k can be bounded 
by (const.)£^|Q;^,fe'|M^^'^~*^'\ The counterterms are fixed in such a way that q;^_oo = 0, 
i.e., a,^- EL-oo(Ci + Fi^^lly' by using the fact that < {const. )elM'^^ and 

< {const. )el\a^^k'\M^^^-^'\ we find that \a^,h\ < (const.)£^*MW2)'». This dimensional 
estimate on a^^h easily implies the desired estimate on it!^(k, p) stated in (4.8) and we will 
not belabor the details here. 



Lowest order computation of a ir^T,T^o+ r^-^A^-^ ^,(2) _ aoi,i,(2) 
is the one-loop contribution to Moreover, (3'^'k^^' = for all k < —1 

neglecting higher order terms, we find ajf^ = ~P'^',o^'^\ ^^^^ is (see Fig. 7): 



Efc<oCfc'^"^ where f3^^^ 



a,l,(2) 

Therefore, 




FIG. 7: Lowest order contribution to ao, ai- 



(2) . -2 

,(0) 



dk 

(2^ 



(2) _ «7l -2 



e 



7.9,„ <7^°)(k + p)C^(k,p)<7(°nk) 



7.9,J^(°)(k + p)C(k,p)^W(k)' 



p=0 



v -'"J (27r)3 
•7.w(°nk). 

After a straightforward computation, using the fact that 
a,J^W(k + p)C(k,p)^W(k)]^^„ = 



Jp=0 



= ^ [ - ^7/.Xo(k) (1 - Xo(k)) + i^df^xoi^)] ^ , 
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(D.9) 



(D.IO) 



where Y = ^o7o + vk-j and (70, 71, 72) = (70, vji, we finally get (3.26)-(3.27). 
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